1 



Prog. Theor. Phys. Vol. ***, No. * May 2011, Letters 

Effects of a New Triple-a Reaction on X-ray Bursts of a Helium 

Accreting Neutron Star 

Yasuhide Matsuo 1 ' *), Hideyuki Tsujimoto 1 , Tsuneo Noda 1 , Motoaki 
Saruwatari 1 , Masaomi Ono 1 , Masa-aki Hashimoto 1 '**) and Masayuki Y. 

Fujimoto 2 

1 Department of Physics, Kyushu University, Fukuoka 812-8581, Japan 
2 Department of Physics, Hokkaido University, 
Sapporo 060-0810, Japan 

o ' 

1 The effects of a new triple-a reaction rate (OKK rate) on the helium flash of a helium 

(— i . accreting neutron star in a binary system have been investigated. Since the ignition points 

determine the properties of a thermonuclear flash of type I X-ray bursts, we examine the 
cases of different accretion rates, dM/dt (M), of helium from 3 x 10 -10 Mo yr _1 to 3 x 
lCP 8 A'/0 yr" 1 , which could cover the observed accretion rates. We find that for the cases of 
low accretion rates, nuclear burnings are ignited at the helium layers of rather low densities. 
As a consequence, helium deflagration would be triggered for all cases of lower accretion 
' rate than M ~ 3 x 1O~ 8 M0 yr _1 . We find that OKK rate is consistent with the available 

observations of the X-ray bursts on the helium accreting neutron star. We advocate that the 
OKK rate is better than the previous rate for the astrophysical phenomena of X-ray burst 
due to helium accretion. 
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2 | 1- Introduction 

■ A new challenge has been given for the triple-a (3a) reaction rate, which has 
been calculated by Ogata et al.® and found to be very large compared with the 
previous rate used so farP^™ As a consequence, the new rate results in the he- 
lium ( 4 He) ignition in the lower density/temperature on the stellar evolution of low-, 
intermediate-, and high-mass stars,®'^ 1 accreting white dwarfs p^ 1 ^ and accreting 
neutron stars.^'^'^ Therefore, it is urgent to clarify quantitatively as possible 

■ as how the new rate affects the above astrophysical phenomena, because the rate 
plays the most fundamental role among the nuclear burning in heavenly bodies and 
could do some role in the early universe, where any terrestrial experiments for the 
3a reaction are very difficult. In the present paper, we investigate the effects of a 
newly calculated 3a reaction rate (OKK rate^ on the helium flashes that occur at 
the bottom layers inside the accreting envelope of a neutron star. We can use the 
ignition curves to find roughly when the nuclear ignition occurs. We note that X-ray 
bursts which have been mainly studied so far are limited to the combined burning 
of hydrogen and heliumP and observational features such as light curves and burst 
energy have been qualitatively explained well. Fujimoto et alP^ have succeeded in 
simulating the X-ray bursts by solving the whole structure equations and clarified 
firstly the importance of the heat flow from the bottom of the accumulated layer. 
The application has been done for several X-ray burst observations of combined H 
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and He burnings ES^SE} However, there remained detailed comparison between ob- 
servations and calculations related to successful bursts, very low accretion rates, and 
or superbursts (e.g.. [T]). |2"3"]) ). 

In the mean time, Fynbo et alP^ (we call the reference as Fynbo) revised the 3a 
rate of NACREP based on new experiments at high temperature of T > 10 9 K and 
with an artificial extrapolation to the very low temperature region toward T ~ 10 7 K. 
However, in the present case of accreting neutron stars we can regard the differences 
between NACRE and Fynbo as unimportant , because the OKK rate is much larger 
for the temperature less than 10 8 K compared to the difference between NACRE 
and Fynbo. 

It is noted that although nuclear burning strongly depends on the temperature, 
the density becomes very important at high densities of p > 10 6 g cm -3 and low 
temperatures of T < 10 8 K due to the screening effects.^ 1 

The new rate has been found to affect significantly the evolution of low-mass 
and intermediate stars 0® where the evolutions from the zero- age main sequence 
through the core He flash/burning from 1 M@ to 10 Mq have been investigated. From 
the HR diagram obtained with the OKK rate, the results disagrees in detail with the 
observations. If the OKK rate in the temperature range of 10 8 K < T < 2 x 10 8 K 
is exactly correct, we must invoke some new physical processes such as rotational 
mixing and/or convection mechanism. 

On the other hand, Saruwatari and Hashimoto^ 1 have shown the important 
difference between the ignition points calculated by the two rates for the helium 
accreting carbon-oxygen white dwarf. They have concluded that for all the accretion 
rates, nuclear fuel ignites at the accumulated layers of helium and a scenario of Type 
la supernovae changes for low accretion rates of M < Ax 10 -8 Mq yr _1 . Related to 
this study, since the progenitor of Type la supernovae have been constructed with 
use of previous reaction rate, discussions concerning the origin of the supernovae that 
originate from the binary between a white dwarf and red giant could be changed.^ 1 

Furthermore, simple helium ignition model on an accreting neutron star has 
been investigated.^ 1 It is concluded that the OKK rate does not fit the observations 
of X-ray bursts due to the pure helium accretion but previous rate agrees with the 
observations. However, their model is based on a model of one variable of column 
density and only steady state has been followed until the ignition. It should be 
noted that the whole structure from the center to the surface of a neutron star 
must be solved to discuss the thermal evolution of accreting neutron stars. This has 
been stressed by Fujimoto et alP^ 1 Therefore, it is desirable to include the whole 
structure and follow the evolution of accreting neutron stars with physical processes 
included,^ if we try to elucidate the effects of the OKK rate on the evolution and 
compare the results with the available observations. We note that for the combined 
burning of H and He, no difference between Fynbo and OKK appears since the 
burst is triggered at higher temperature region of T > 5 x 10 8 KP In this letter, 
we present the results of the evolutionary calculation beyond the helium ignition 
of a helium accreting neutron star. We will show the important role of both the 
3-a reaction rates in the low temperature region and the properties of the crustal 
heatingPSJ'E) 
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2. Evolution toward the explosive helium burning on the accreting neutron star 

Accreting neutron stars are considered to be the origin of the type I X-ray 
bursts.^ Accreting materials are usually hydrogen and/or helium. Since the hy- 
drogen is converted to helium through steady hydrogen burning, helium is gradually 
accumulated on the neutron star and the deep layers become hot and dense. Helium 
flash triggered in the region composed of degenerate electrons could develop to the 
dynamical stage, depending on the accretion rate M.^ 1 The properties of ignition 
are determined from the thermal structure around the bottom of accreting layers. It 
has been found that unstable helium burning after the exhaustion of hydrogen results 
if the helium accretes in the range of 2 x 10 _10 M Q yr _1 < M < 1O -9 M yr _1 P3 

We have computed the evolution of an accreting neutron star by the Henyey- 
type implicit-explicit method-^ 1 during the hydrostatic evolutionary stage of the 
accretion of helium, where the full set of general relativistic equations of spherically 
symmetric stars is based on the formulation by ThorneP^ 1 Physical inputs are the 
same as Ref I19p except for the network of nuclear burning and additional heating 
from the crust. We note that under the assumption of the spherically symmetry, the 
basic four equations of neutron star evolution becomes the same type of those of the 
ordinary stellar evolution except for the effects of the general relativity. 

Initial models are constructed by accreting helium on a neutron star of which 
the gravitational mass and the radius of the neutron star are M = 1.3 M Q and R = 
8.1 km, respectively. We can obtain the initial models for each constant accretion 
rate shown in Table HI after the steady state has been achieved as in Ref. fl~9j) : the 
integrated energy flux of the non-homologous part for the gravitational energy release 
becomes negligible and therefore the thermal structure does not almost change. The 
masses have been gradually increased during the accretion toward the steady state 
and the amounts are less than 1 % of the total mass when the steady state has been 
achieved. 

To follow the nuclear burning, the alpha-network has been implemented to obtain 
the nuclear energy generation rateP3 the network consists of 4 He, 12 C, 16 0, 20 Ne, 
24 Mg, 28 Si, 32 S, 40 Ca, 44 Ti, 48 Cr, 52 Fe, and 56 Ni. We include and examine the 
two reaction rates of Fynbo and OKK for the 3a reaction, and the other reaction 
chain is selected from the faster reaction of either (a, 7) or (a,p) followed by the 
(p, 7) reactions. After the calculations of stellar structure have been converged, 
this network is operated to obtain the nuclear energy generation rate to calculate 
the next stellar structure as the same method in Ref. I26|) . The nuclear energy 
generation rate of the 3a reaction is crucial to determine the ignition condition and 
would be proportional to T g -3 exp(— 44.0/Tg^P3 for the high temperature region of 
Tg > 0.1. However, the temperature dependence cannot be written by a simple 
analytical formula, that is, 

esaocp 2 ^ f(3a)(T), (1) 

where Y is the helium mass fraction and / is the screening factor which is included 
by the same method in Ref. [5]). As a function of T, (3a) (T) was obtained from 
the average with the Maxwell-Boltzmann distribution concerning the product of the 
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Fig. 1. Evolutionary track of the temperature against the density at the bottom of the burning 
layer in M = 3x 1O _1O M0 yr _1 for the 3a reaction rate of OKK (left panel) and Fynbo (right 
panel). The ignition and deflagration curves are defined by e n = e ra< j and Eq.lJSJ, respectively. 



cross section and the relative velocity.® 

The importance of a crustal heating during the thermal evolution of neutron 
stars has been pointed outP^ The heating rate is tabulated in Ref. |28|) and the 
heating rate is 

Qi = 6.03M_i ^j|^ 10 33 erg s" 1 , (2) 

where M_io is accretion rate in units of 10 -10 Mq yr -1 and qi is deposited heat in 
MeV/nucleon whose representative value within the present accretion rate in this 
paper would be evaluated to be 0.012 MeV/nucleon from the produced average 
heating rate: 

<Qi>= / Qidm/A M, AM = / dm, (3) 

J crust J crust 

where dm = 4nr 2 pdr and the integral covers the crust of the density range 10 11 — 
10 13 g cm -3 and A M ~ 1.2 x 10~ 4 M with the thickness of 1 km. We can easily 
include it just as the neutrino energy loss rates and/or nuclear generation rates, 
that is, the rate of the crustal heating is added in addition to the nuclear energy 
generation rate, where we write the energy equation in Newtonian approximation 
for simplicity: 



^--e* -e -t(— 

dM r ~ £n £v \dt 



where L r is the energy flow, e* includes the nuclear energy generation rate, crustal 
heating rate, e v is the neutrino energy loss rate, s is the specific entropy. 

We have selected six cases of accretion for two reaction rates (see Table |T]) , where 
the hight of the accreted layer is at most some 10 m during the shell flash. We have 
shown the evolutionary tracks in Fig. Q] for the OKK rate (left panel) and the Fynbo 
rate (right panel) with the solid lines. The tracks correspond to the bottom of the 
accretion layers, where the maximum temperature results and shell flash should be 
triggered in the neighborhood of the region. 
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Let us define the ignition and deflagration and/or detonation curves to show 
clearly the points of the helium ignition and the development of the shell flash. The 
density and the temperature in the below equations (4)-(8) must correspond to the 
bottom of the burning layer and/or around that of the accreted layer. The energy 
conservation law can be written as follows; 

CP — = En ~ £rad, (4 J 

where cp is the specific heat at the constant pressure and e n is the nuclear generation 
rate of the 3a reaction. We approximate the radiative loss in the simple form, 



4acT 4 



where k is the opacity, a is the column density, a is the radiation constant, and c 
is the speed of lightP^ 1 The column density can be obtained from the pressure and 
the gravitational acceleration for a neutron star model (see Eq. ([9])). The ignition 
curve is defined on the plane of (p,T) to satisfy the equality e n = e ra d- The ignition 
curves are drawn in Fig. [1] for the OKK rate and the Fynbo rate, where we can infer 
that helium ignition should be triggered at the much lower density. 
The dynamical timescale Td yn is defined as 

Tdyn = (6) 

where G is the gravitational constant. On the other hand, the time scale r n of the 
increase of the temperature due to the nuclear burning is defined to be 

Tn = ——■ (7) 

En 

The deflagration curve has been obtained from the below condition^ 

Tn = Td y n- (8) 

If t ti < 7d yn , deflagration wave should be originated, which could develop and burn 
out the previously accumulated layers (hereafter we call the wave the deflagration 
one). The deflagration curves are also drawn on the plane of (p,T) in Fig. [TJ We 
note that the curves depend significantly on the helium mass fraction designated by 
Y. 

Figure 1 shows the evolutionary track on the plane of (p, T) toward the helium 
ignition for the bottom of the burning layer (around the bottom of the accreted layer) 
and subsequent helium flash beyond the deflagration curves, where the deflagration 
curves of two cases of helium mass fractions (Y = 0.1 and 0.5) are drawn with use of 
the OKK rate and Fynbo rate. We can find that the helium ignition occurs at lower 
density regions by almost two orders of magnitude if the OKK rate is adopted. We 
note that for M < 3 x 10 -8 M Q yr _1 , the helium flash develops to the dynamical 
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Table I. Energy releases -Eburst per a burst for a fixed accretion rate are given for two cases of OKK 
and Fynbo. The ignition pressure for each M , Pi gn , corresponds to the maximum temperature 
layer inside the accretion layers whose location is around the bottom of the accumulated ones. 
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stage. Therefore, our hydrostatic evolution code cannot follow the dynamical stage 
which would lead to the deflagration and/or detonation. 

3. Comparison with the observations 

It would be useful to judge which 3a rate is more consistent with available obser- 
vations. In general, it is very difficult to follow in the hydrodynamical calculations 
the helium flash until the fuel of helium depletes as described in the previous section. 
However, using the results of evolutionary calculations obtained in the previous sec- 
tion, we can compare qualitatively the theoretical calculations with the observations. 
We choose the layer of the helium ignition from which the pressure, the density, the 
temperature, and the composition of helium are obtained. We define the density 
Pign and the pressure Pi gn for the layer when the helium ignition occurs. We can 
get 

Agn from the crossing point between the evolutionary track of the bottom of the 
accreting layer and the ignition curve. As a consequence, Pi gn is obtained from the 
equation of state. Since the accumulated layers are well approximated by the plane 
parallel configuration, we can adopt the approximation to estimate the accumulated 
mass with a enough approximation; It is written in the spirit of the plane parallel 
model as follows PB >EQ ,E3 

-Pign = 9s°~igm (9) 

where g s is the gravitational acceleration and <7i gn is the the column density at the 
ignition point. The energy release -Eburst can be obtained under the assumption that 
the accumulated layers on the ignition layer are completely burnt out: 

£ b urst = AM ign Q nuc /{I + z s ), (10) 
= 47r^ 2 a ign Q nuc /(l + z s ). (11) 
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Fig. 2. Energy release of -Eburst against the accretion rate M for two cases of OKK and Fynbo. 
Observational values from different astrophysical objects which have been identified to cause 
type I X-ray bursts due to the pure helium accretion, are plotted as was done in Ref. I15|>. 



Here, AM[ gn is the total rest mass of the accreted layers, the gravitational redshift, 
z s = (1 — 2GM t /c 2 R)~ 1 / 2 — 1 ~ 0.38, and Q nuc is the nuclear energy release per 
nucleoli ((5 nU c = 1.6 MeV/nucleon), where all fuel of helium are assumed to be 
burnt into iron. For our model of the neutron star, since g s = GM(1 + z s )/R 2 , we 
get log # s = 14.56. 

In Table [H we show Pi gn and -Eburst in each accretion rate for the cases of OKK 
and Fynbo. In Fig. we show Eburst as a function of accretion rates M, where are 
shown the two observations of type I X-ray bursts, 4U1820-30 (labeled by 4U1820- 
30), SLX 1737-282, SLX 1735-269, and 2S 0918-549 (labeled by intermediate long 
bursts). Their accreted matter from companions has been known to be pure he- 
l[ um ^XMbEB£MX3BMb As is written in Ref. [15]), the uncertainties of accretion 
rates and burst energies are indicated by boxes. For these three observations, we 
can see from Fig. [21 results by OKK are consistent with the observations if we choose 
the accretion rate of M < 3 x 1O _1O M0 yr -1 . 

4. Discussion and conclusions 

The ignition densities that determine the triggering mechanism on helium ac- 
creting neutron stars will be changed significantly if we adopt the new 3a reaction 
rate for the low accretion rate of M < 10 _8 Mq yr . Until the OKK rate has been 
presented, the method by Nomoto^ 1 has been used of which a simple extrapolation 
of the Breit-Wigner type function to the low-temperature side is adopted, where it 
has been advocated that the non-resonant 3a reaction is crucial in determining a 
helium ignition of compact stars for the low accretion rate. However, the microscopic 
calculation for the three body problem has been found to be crucial in evaluating 
the 3a reaction rate. In the present paper, we have examined the effects of the OKK 
rate on the very low temperature site of astrophysics. Although our results seem to 
contradict with those of Peng and OttJ^ we may ascribe it to the following. First, 
we have adopted the rate of the crustal heating from Ref. [28|) . Second, they have 
used the mass of 1.4 Mq and the radius of 10 km for a neutron star. Third, while 
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M (M yr- 1 ) 



Fig. 3. Energy flux in units of erg cm 2 s 1 from the crust against accretion rates with use of the 
tabled (solid line) and the two boxes are taken from Ref. I15jl . 

we have followed the evolution of an accreting neutron star from the center to the 
surface of the star, they have solved only steady-state equations above the crust as 
a function of a single dependent variable of the column density. As a consequence, 
it has been written that the outward energy release due to the crustal heating was 
taken to be a free parameter as a boundary condition, which would be inappropriate 
in the calculations of accreting neutron stars. 

The most important difference comes from the treatment of the crustal heat- 
ing. In the present case, the energy flux (erg cm~ 2 s~ 1 ) from the crust is shown in 
Fig. [3l The average energy per unit mass from the crust can be evaluated to be 
~ 0.01 MeV/nucleon as shown in the previous section. On the contrary, the rate 
they have adopted could be nearly 0.1 MeV/nucleon These differences reflect the 
calculations of the heat flux as shown in Fig. O where the solid line represents heat 
flux in our calculations and the two boxes are taken from Ref. I15p . It is clear that the 
heat flux from the crust in the present study is 1-2 orders of magnitude lower than 
the case in Ref [T5|) . As a consequence, our initial models just before the accretion 
starts might have lower temperatures compared to those used in Ref. I15p . The most 
significant difference would come from the initial model and therefore, the location 
of the ignition point on the (p, T) plane. While the bottom of the accreted layers in 
our initial model for M = 3 x 10 _10 Mq yr _1 has the temperature of T = 1.6 x 10 K 
(see Fig. [TJ), we can infer for the case of Ref. [T5|) the temperature could be rather 
high. From the column density obtained from the box of the observation of interme- 
diate long bursts (see Fig. 2 in their paper), we can get the ignition temperature of 
T = (1 — 2) x 10 8 K for the Fynbo rate. If the initial temperature of their calculation 
take nearly the same temperature, which is the common feature of the accreting 
compact objects before the beginning of the shell flashes ignition density is found 
to be in the neighborhood of that for the OKK rate. Contrary to our calculations 
which are the results of stellar evolution and include all the necessary physical pro- 
cesses, their calculations depend essentially on the numerically convenient boundary 
condition at the assumed surface of the crust. 
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We have shown that OKK rate must change the scenario of Type la supernovae 
and type I X-ray bursts for the low accretion rates, where the temperature at the 
ignition is less than 10 8 K. Since the rate was calculated under the situation of three a 
particles in vacuum, application of the OKK rate is doubtful against the high density. 
As seen from Fig. 3 in Ref . H]) , the discretized continuum wave functions of the two 
a particle system extend to 5 x 10 3 fm. Therefore, in the medium of stellar plasma, 
effects more than three a particles cannot be negligible for p > 10 6 g cm" 3 P For 
example, the reaction rate is to be calculated with the inclusion of an appropriated 
screening potential. If the effects of many body interactions are included, the OKK 
rate could be decreased by orders of magnitude. We may constrain the 3a reaction 
rate from the astrophysical observations of type I X-ray bursts. 

As example of constraint the OKK rate, if we adopt the upper values of the 
observational box of .E burst in Fig. [2]for M = 3 x 1CT 1O M yr 1 , we obtain log P\gn ~ 
24.9 from ([9|)- ()lip . Since the evolutionary paths are almost the same for the two 3a 
rates until the ignition points (see Figs. [1]), the corresponding ignition density with 
the evolutionary paths of OKK becomes in the range of \ogp lsa = 7.7 — 7.8 for the 
OKK rate reduced by a factor of 10 2 ~ 3 . We note that the properties of shell flash 
depend to some extent on the structure of a neutron star (mass and radius); Since 
^burst is proportional to R 4 /M, a rather hard EOS (R = 12 km and M = 2.0 M ) 
gives three times larger value compared to our case, which would be preferable for the 
OKK rate. On the other hand, it depends on the non-standard cooling processes: 
if we use the non-standard cooling processes such as a pion-condensationj^ our 
estimates would be changed significantly. 

In conclusion, we have found that the new rate affects the helium ignition in ac- 
creting neutron stars, where for lower accretion rates, helium burns at lower densities 
and temperatures, which should change the epoch of the formation of a helium de- 
flagration wave and the modeling of type I X-ray burstsP In particular, we consider 
that the crustal heating affects the ignition properties significantly. Furthermore, the 
mechanism behind superbursts should be studied again using the new rate, because 
the amount of 12 C could play a crucial role in inducing the superbursts.^ 
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